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The properties of the two-flavored Gross-Neveu model in the (l+l)-dimensional x spacetime 
with compactifled space coordinate are investigated in the presence of the isospin chemical potential 
Hi. The consideration is performed in the limit Nc oo, i.e. in the case with infinite number 
of colored quarks. It is shown that at L = cxd (L is the length of the circumference S^) the pion 
condensation phase is realized for arbitrary small nonzero /ij. At finite values of L, the phase 
portraits of the model in terms of parameters v ^ Hi and A ~ 1/L are obtained both for periodic 
^•ZJ ' and antiperiodic boundary conditions of the quark field. It turns out that in the plane (A, v) there 

, is a strip < A < Ac which lies as a whole inside the pion condensed phase. In this phase the pion 

condensation gap is an oscillating function vs both A (at fixed u) and u (at fixed A). 
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INTRODUCTION 



' It is well-known that QCD is a fundamental theory of strong interactions both in the vacuum and in hot and/or 
dense baryonic matter. However, it can be successfully used only in the region of high energies, temperatures and 
^ • densities (or chemical potentials), where a weak-coupling expansion is applicable. Away from this region, different 
(]J ' nonperturbative methods or effective theories such as chiral effective Lagrangians as well as Nambu - Jona-Lasinio 
. type models (see, e.g., the papers and references therein), are usually employed for the consideration of light 

meson physics, phase transitions in dense quark matter, etc. In particular, motivated by the fact that in heavy-ion 
, collisions and compact stars the hadronic matter is isotopically asymmetric, different QCD-like effective models were 
^ ' studied at nonzero isospin chemical potential /x/ @, @, 0, S Q . There the charged pion condensation phenomenon, 
, which is generated if /i/ is greater than the pion mass m^, was also considered. 

■ In all above mentioned papers the effective models are i) field theories in usual (3-|-l)-dimensional spacetime, and 
ii) they are employed for the description of QCD at rather low energies and densities. At the same time there 
is another class of theories that can be used as a laboratory for a qualitative consideration of QCD at arbitrary 
' energies. These are the so-called Gross-Neveu (GN) type models, i.e. two-dimensional quantum field theories with 
^5 , four-fermion interactions fldl . [Til . [T^ . Renormalizability, asymptotic freedom as well as the spontaneous breaking 
of chiral symmetry (in the vacuum) are the most fundamental features that are inherent both for QCD and all GN 
type models. In addition, the GN phase portrait in terms of baryon chemical potential I-Ib vs temperature resembles 
qualitatively to a great extent the QCD phase diagram [l^ 0, [H, [H, [13] ■ Due to their relative simplicity in the 
leading order of a large A'c-expansion (Nc is a number of colored quarks), it is very convenient to use GN models for 
considering such a phenomenon of dense QCD as color superconductivity p^[l8| and to elaborate new nonperturbative 



methods of quantum field theory [19, 20, 2lJ. Moreover, the influence of the space compactification on chiral symmetry 



22 1 and in dense baryon matter (/is 7^ 0) [23] was studied in terms of GN 
24 



breaking both in the vacuum {^b — 0) 
models (see also the appropriate papers \24. l25l. l26j|l. 

Before investigating different physical effects relevant to a real (3+l)-dimcnsional world in the framework of two- 
dimensional GN models, let us recall that there is a no-go theorem forbidding the spontaneous breaking of continuous 
symmetries in two dimensions (27| . However, at present time it is well understood (see, e.g., the discussion in 
[ill , [isl m, [13) that in the limit Nc 00 this no-go theorem does not apply. This makes it possible to study 
symmetry breaking effects in terms of GN models as well, but only in the leading order of the l/7Vc-expansion, where 
most low dimensional theories are exactly solvable. In this sense, for Nc — s- 00, low dimensional quark models are 
physically more tractable and appealing than at finite Nc- 

In the present paper the pion condensation phenomenon is investigated in the framework of the two-dimensional 
GN model with two massless quark flavors. In particular, we shall study the influence of the finiteness of the system 
size on this phenomenon. So our consideration is performed in a spacetime with non-trivial topology, i.e. on the 
X manifold with compactifled space coordinate, and the GN model is extended by an isospin chemical potential 
fii (for simplicity, we put iib = 0). Obviously, the latter issue is motivated by the physics of compact stars, where 
pion condensation might be realized as a consequence of the isotopic asymmetry of baryon matter. Since all the 
calculations are carried out on the basis of the leading order of l/7Vc-expansion (i.e. in the case Nc 00) we 
expect that all conclusions concerning the pion condensation phenomenon, caused by a spontaneous breaking of the 
continuous isospin symmetry, remain qualitatively valid for real QCD. 
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II. THE CASE OF x SPACETIME 
The model and its thermodynamic potential 



We consider a two-dimensional model which describes dense quark matter with two massless quark flavors {u and 
d quarks). Its Lagrangian has the form 



L 



i^-o , Ms , A*/ 
7 lOi' + —7 + Y'r37 



G 



(1) 



where the quark field q{x) = qia{x) is a flavor doublet (i — 1,2 ov i — d) and color Nc-^\et {a — 1, Nc) as well as 
a two-component Dirac spinor (the summation in ([T]) over flavor, color, and spinor indices is implied); (fc — 1, 2, 3) 
are Pauli matrices; the baryon chemical potential ^ib in ([T|) is responsible for the non-zero baryon density of quark 
matter, whereas the isospin chemical potential /i/ is switched on in order to study properties of quark matter at 
nonzero isospin densities (in this case the densities of u and d quarks are different). Evidently, the model ([1]) is a 
generalization of the two-dimensional Gross-Neveu model with a single massless quark color TV^-plet to the case 
of two quark flavors and additional chemical potentials. As a result, we have in the case under consideration a more 
complicated chiral symmetry group. Indeed, at /z/ = apart from the global color SU(A^c) symmetry, the Lagrangian 
([1]) is invariant under transformations from the chiral SUl{2) x SUii(2) group. However, at /i/ 7^ this symmetry 
is reduced to Ui^l{^) x Ui^fi{l), where I3 = T3/2 is the third component of the isospin operator (here and above 
the subscripts L,R mean that the corresponding group acts only on the left, right handed spinors, respectively). 
Evidently, this symmetry can also be presented as Ui^{l) x [7^/3(1), where Ui^{l) is the isospin subgroup and UAisi^) 
is the axial isospin subgroup. Quarks are transformed under these subgroups as q ^ exp(iQ;T3)g and q — > exp{ia'~f^T3)q, 
respectively. ^ 

The linearized version of the Lagrangian which contains composite bosonic fields a{x) and iTaix) (a — 1,2,3), 
has the following form (in what follows, we use the notation /i = /iB/3 for the quark chemical potential): 



L 



7'^ia. + M7° + ^r37" 



cr - 17 TTaTa 



4G 



aa + TTaTTn 



From the Lagrangian ([2]) one gets the equations for the bosonic fields 



G 

a{x) ^ -2 — {qq); TTa{x) 



G 

-2^(9i7'^a'7). 



(2) 



(3) 



Obviously, the Lagrangian Q is equivalent to the Lagrangian ([T]) when using the equations ([3]). Furthermore, it is 
clear from ^ and footnote [T] that the bosonic fields transform under the isospin 1/73(1) and axial isospin [74/3(1) 
subgroups in the following manner: 



[7/3(1): a a; tt^ tt^; tti ^ cos(2a)7ri + sin(2a)7r2; 7r2 
[7^1/3(1): TTi ^ TTi; TT2 ^ TT2] (T cos(2q!)(t + sin(2a)7r3 ; 



cos(2a)7r2 — sin(2Q)7ri, 
cos(2a)7r3 — sin(2a)fT. 



(4) 



Starting from the theory one obtains in the leading order of the large A'^c-expansion (i.e. in the one-fermion loop 
approximation) the following path integral expression for the effective action iScff (cr, tTq) of the bosonic (7{x) and na{x) 
fields: 



where 



exp(i5off(cr,7ra)) = N' I [dq] [dg] cxp ( i 



j d^x 




4G 





Ld^x 



(5) 



N' is a normalization constant. The quark contribution to the effective action, i.e. the term ScS in ([5]), is given by: 



exp(i5cff) = N' J[dq][dq]exp(^i J 



qDq 



d'^x) ^ [BetD] 



(6) 



In ^ we have used the notation D = Z) x Ic, where Ic is the unit operator in the Ac-dimensional color space and 

D = 7'^i5, -I- M7" + Y^3l" '^l^^aTa 



(7) 



Recall for the following that exp(iar3) = cos a + ira sin a, exp(iQ7^T3) = cos a + i7^T3 sin a. 
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is the Dirac operator, which acts in the flavor-, spinor- as well as coordinate spaces only. Using the general formula 
DetD = expTrlnZ?, one obtains for the effective action the following expression 



5off((T, TTa) 



Nc / (fx 



4G 



(8) 



where the Tr-operation stands for the trace in spinor- (s), flavor- (/) as well as two-dimensional coordinate- (x) spaces, 
respectively. Using ([5]), we obtain the thermodynamic potential (TDP) fi^^^^ (cr, tTq) of the system: 



cr, TTa— const 
72 



AG 



. Tisfx In D 



Jdh 



AG 



iTr 



In ^ + /i7 + YT-37 - Cr - 17 TTaTa 



(9) 



where the a- and ttq fields are now x-independent quantities, and in the round brackets of ([9|) just the momentum 
space representation, D, of the Dirac operator D appears. Evidently, Tr^/ In I? =Y^^ Ine^, where the summation over 
all four eigenvalues ti of the 4x4 matrix D is implied and 



ei,2,3,4 = -cr ± W (po + /i)^ -pI~'TtI + (m//2)2 ± A^/V (po + m)^ - -^l 



(10) 



Hence, 



4G 

^2 



(2^ 



2 



where 



PI 



3 2 



(11) 



(12) 



The TDP f2p,^j (ct, TTa) is symmetric under the transformations /i —t — /i and/or /i/ —t —fii. Hence, it is sufficient to 
consider only the region /i > 0,/^/ > 0. In this case, one can integrate in (|lip over pq with the help of the formula 



dpo 
2ti 



In 



{p^ + af-b"] = '-[\a-b\ + \a + h\] 



(13) 



(which is valid up to an infinite constant independent of quantities a, h) and obtain: 



• {cr, Ha) 



AG 
AG 



\e-+^l\] 



— [e+ + e- + {^i- e+)e{^i - e+) + e-)e{^i - £_)|. 



r 



(14) 



(To get the second line in we used the relations |£± + /i| = £± + /i and d{x) + 0{—x) = 1.) In what follows we 
are going to investigate the /i, /i/-dependence of the global minimum point of the function (cr, tt^) vs cr, tTo. To 
simplify the task, let us note that both the quasiparticle energies (|12p and hence the TDP (I14p depend effectively 
only on the two combinations cr^ + 7r| and ttJ + 7r| of the bosonic fields, which are invariants with respect to the 
L//3(l) X UAi^i^) group, as is easily seen from (|4]). In this case, without loss of generality, one can put 7r2 = tts = in 
([T^ . and study the TDP as a function of only two variables, M = a and A = tti. Then the global minimum point of 
the TDP 17^,^, (M, A), 



AG 



r ^i^^ + + (a* - E+)e{^i - E+) + (^ - E-^)e{^i -e-^)], 



(15) 



is the solution of the system of gap equations 



9r!^,^,(M,A) ^ M 
dM ~ 2G 
dn^,^,{M,A) _ A 



M 



d^r0{E+- fi)E+ , e{E^-fi)E- 
2TrEV 



El 



E-, 



dA 



2G 



dp^ce{El^ e{E^-_jj^ 



2ti I 



El 



E-, 



(16) 
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where = ^/{E^Y + A2, E^ = E±^, !md E = ^Jpl + JVP. Evidently, the coordinates M and A of the global 
minimum point of the TDP (|15p supply us with two order parameters (gaps), which are proportional to the ground 
state expectation values of the form {qq) and (qi'y^Tiq), respectively. If the gap M is nonzero, then in the ground 
state of the model the axial isospin symmetry UAi-ji^) (at /i/ ^ 0) is spontaneously broken down. Moreover, if the 
gap A ^ 0, then in the ground state, corresponding to the phase with charged pion condensation, the isospin Ui^(l) 
symmetry is spontaneously broken down. 



B. Pion condensation: the case of jj, — 0, jii 



Since at uj = 0, /i 7^ the phase structure of different GN models was reasonably well studied both in two dimensions 
[isl . [isl . [ig . and in three dimensions [l^ (in the last case the four-fermion theories are also renormalizable) , in 
this subsection we shall study for simplicity the model ([1]) only at zero quark chemical potential, i.e. at /i = 0, but 
fii ^ 0. The corresponding TDP will be denoted as fl^^^M, A) and can be obtained from pS]): 



Voip)~ J^^^{e++E^-2^p^+pI], (17) 



where p — V A'P + A^ and Vo{p) is the TDP of the system in the vacuum, i.e. at pi — 0. In the vacuum the TDP is 
usually called effective potential: 

It is easily seen that both the TDP (fT7|) and the effective potential (HH]) are formally ultraviolet (UV) divergent 
quantities. So, a few words are needed about the renormalization procedure of the initial model. It is well known 
that all four-fermion theories of the type (H]) are renormalizable in two dimensional spacetime Moreover, in the 
leading order of the large Nc expansion only the coupling constant should be renormalized in order to obtain finite 
(renormalized) expressions for different quantities (see, e.g., [HI). It means that the bare coupling constant G of the 
model (HI) depends on the cutoff parameter A, G = G(A), in such a way that all UV divergences, arising from loop 
integrations when A ^ 00, are compensated by corresponding terms of G(A). As a consequence, in the limit A ^ 00 
one must necessarily obtain finite expressions for physical quantities. Of course, different renormalization procedures 
result in different expressions for the bare coupling constant G(A). However, physical consequences of the theory do 
not depend on the concrete renormalization scheme. Taking this last remark into account, let us next discuss how to 
obtain finite renormalized expression for the TDP (|17p . We see here two ways. On the one hand, one could find an 
expression for the bare coupling constant G such that the UV divergence, arising from the integral in the first line 
of p?)) . would be compensated by the term with G. Evidently, in this case G depends both on the cutoff A and pj. 
However, we find it more convenient to consider the second way. In this case, one should first of all note that the 
integral in the second line of pT|) is a convergent quantity, and the whole UV divergence is located in the effective 
potential Vo(p). Hence, there is a possibility to remove UV divergences using a bare coupling constant which does 
not depend on pi. Namely, let us choose 

where Mq is the dynamical mass of quarks in the vacuum (for more details, see Appendix E]) . Then, substituting 
(fro)) into fTS]) and restricting there the range of integration by using the cutoff parameter A, it is possible to obtain 
for A >> Ma the expression (moreover, we omit an inessential infinite constant independent of p): 



Voip) = f 



In 



(20) 



Since AIq might be considered as a free model parameter, it follows from (ITni) - (|20p that the renormalization procedure 
of the GN model is accompanied by the dimensional transmutation phenomenon. Indeed, in the initial unrcnormalized 
expressions both for 51^^(1/, A) and Vo{p) (see (IT7)) and (fTS]) . respectively) the dimensionless coupling constant G 
is present, whereas after renormalization the effective potential (|20[) is characterized by a dimensional free model 
parameter Mq. 

Due to the relation ([20]) . one can show that the gap equations for the renormalized TDP (fT7|) might have no more 
than three different solutions. Two of them, (M = 0, A = 0) and (M = 0, A = Mq), are present at arbitrary values 
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of III > 0, whereas the third one, (M = Mq, A = 0), appears only at /i/ < MqV^- However, for arbitrary /i/ > a 
global minimum point of the TDP r2^j-(M, A) lies at the point (M = 0, A = Mq). This means that in the model ^ 
the isospin symmetry is always broken down and a charged pion condensate which is equal to the quark mass Mq in 
the vacuum, is created if /i/ > 0. 

Since in the vacuum case (/i — 0,fii ~ 0) chiral symmetry is spontaneously broken down in the model ([1]), there 
must exist three massless Nambu-Goldstone bosons which are pions, i.e. = 0. So, we have proved that in the 
framework of the model (H]) the pion condensation phase is realized at fij > rriT^, where m-^ is the pion mass in the 
vacuum. Just the same phase structure is predicted by QCD at /i = 0, /z/ 7^ [5|. In contrast, in the framework of 
(3+l)-dimensional NJL-type models the pion condensation is not allowed for sufficiently high values of the isospin 
chemical potential @, 0, @| • This fact supports the statement made in the Introduction that the N JL approach is only 
valid at rather small energies (chemical potentials). Moreover, we have once more demonstrated that in the leading 
order of the large A^c-expansion the two-dimensional GN models are a quite good theoretical laboratory for qualitative 
QCD investigations. So we are in a position to believe that the results obtained in the next sections are also inherent 
to QCD. 



III. THE CASE OF x SPACETIME AND m/ / 

In the present section we continue the investigation of the charged pion condensation, this time under the influence 
of the finite volume occupied by the system. This is obviously a reasonable task, since all physical effects take place 
in restricted space regions. The consideration of the problem is significantly simplified in the framework of the two- 
dimensional model ll]) at 7^ 0, which is again justified by its similarity to QCD. So we put a system with Lagrangian 
([T|) into a restricted space region of the form Q < x < L (here x is the space coordinate) . It is well known that in this 
case the consideration is equivalent to the investigation of the model in a spacetime with nontrivial topology x S*^ 
and with quantum fields, satisfying some boundary conditions of the form 



q{t,x + L) ^ e"°'q{t,x), 



(21) 



where < a < 2, L is the length of the circumference S^, and the variable x means the path along it. Below, we shall 
use only two values of the parameter a: a = for periodic boundary conditions and a = 1 for the antiperiodic one. 

As a consequence, to obtain the thermodynamic potential Q^^^^^M, A) of the initial system placed in the restricted 
domain < x < L and at /x/ 7^ 0, one must simply replace the integration in p7p by an infinite series, using the rule: 



dp 



-f{Pi) 



-{2n 



0,±1,±2,. 



(22) 



n— — 00 



Moreover, instead of Vo{p) it is necessary to use the effective potential Vl{p) of the model in the vacuum (see Appendix 
E]) . As a result, the TDP ([T7]) will be replaced by the corresponding expression for the spacetime of the form x S^, 



n— — 00 I \ \ 



2„ + a)2 + ^ +A2 



\ 



'M2 + ^(2n + a)2-^ 



A2-2\/p2 + _(2n + a)2 



(23) 



where p = V M2 -|- A2, and the function Vl{p) is defined in (jAlOp . In what follows, it will be convenient to use the 
dimensionless quantities 



A 



LMo 



2Mo' 



M 

Mo' 



(24) 



where AIq is the dynamical quark mass in the vacuum. Moreover, since the phase structure of the model in the two 
particular cases L = 00, pi ^ and L ^ 00, pi = was already considered in section FlI Bl and in the Appendix [Xl 
we will now investigate the phase structure only at A > 0, > 0. 
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FIG. 1: The periodic case: In the left picture the phase portrait of the model is represented in terms of (A, v), where number 
1 denotes the symmetric phase, number 2 - the pion condensed phase. In the right picture the behavior of the gap (5 vs A is 
depicted at i/ — 7.5. 

A. The case of periodic boundary conditions 

In this case a = 0, and in terms of the dimensionless quantities (p4)) the TDP ((23)) can be rewritten in the following 
explicit form 



- 2A I [yjm^ + (2nA)2 + i/) + '^^ + y (v^^^ + ^ t^) " + - 4nA - ^ ' | , (25) 

where 7 — 0.577... is the Euler's constant [29i]. We consider the TDP (P5|) as a function of two variables, m,(5. 
Moreover, A are free parameters there. Since the information about the phase structure of the model in the case 
of the periodic boundary conditions is contained in the global minimum point of the function (|25p vs m, 5, it is first 
of all necessary to study the gap equations and then to investigate the behavior of the global minimum point vs 
parameters A. In particular, it is possible to show that for each fixed point of the plane (A, v) (with z/ > and 
A > 0) the global minimum point of the TDP ([^5]) might be located at two different points only, i) (m = 0, (5 = 0) 
and ii) (m = 0, (5 7^ 0), where the nonzero gap 8 is the solution of the equation 9Daj/(0, (5)/9(52) = 0. The point 
i) corresponds to the Ui^(X) x Uah^X) symmetric phase of the model (without charged pion condensation). On the 
other hand, if the global minimum of the function (|25D is situated at the point ii), then in the ground state of the 
model the isospin symmetry [//^(l) is spontaneously broken down, and the pion condensation takes place. Let us 
denote by Ic the critical curve which separates the region of the (A, v) plane with symmetric phase from the points 
(A, v\ corresponding to the pion condensed phase of the model. Since in each point of the curve /c there is a phase 
transition of the second order from the symmetric phase to the pion condensed one and vice versa, the gap 5 must 
vanish on this curve. So the critical curve Ic is defined by the following equation 

9DA,.(m, (5) 



A °° r A A 11 

9(^2) - - 7 - - - E 1^;^ + - - j - 0. (26) 

To represent the curve Ic in the plane (A, v)^ it is convenient to divide this plane into an infinite set of regions cj^: 

00 

(A,i/)=|Jwfc; Wfe = {(A,!/) : 2A(fc- 1) < < 2Afc}. (27) 
fc=i 

C30 

In accordance with the division (P7)) . the critical curve 1^ can also be presented as a set of pieces, /c = U ^cfc- 

fe=i 

Obviously, each piece l^k of the whole critical curve lr_ lies inside the corresponding fc-th region uj^ and obeys the 
following equation (k > 1) 



A A 1 



''ck 



2nA + V 2nX - 



v n 



- - = 0. (28) 
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„J . . . . , . . . . , . . . . , . . . . , . .J . . . . , . . . . , . . . . , . . . . , . . . . , rl/ 

" 2468 "2468 10 

FIG. 2: The periodic case: The behavior of the gap 5 vs at A = 1 (left picture) and A — 1.7 (right picture). 

For k — 1 the part Id obeys the equation (|26p with omitted absolute value symbols. Performing the summations in 
((26)) or ([28]) . one can find for each piece Ick of the critical curve Ic the following equation (fc > 1) 

Ick-. 21n(4A) + 2^(/c-^)-^.(l-^)+^(^) -0, (29) 

which is valid only at 2A(fc — 1) < < 2Xk. Here ip{x) is the logarithmic derivative of the Euler's r{x) function [2^. 
Before drawing the critical curve Ic, we would like to point out one its peculiarity. Using the well-known property of 
the ipix) function, 7rcot(7ra;) = — x) — 'ijj{x), as well as the periodicity of cot(7rx), the equation pSj) can be reduced 
to the following one: 

Ick ■■ 2 ln(4A) = -2iIj{z) - tt cot(7rz) = F{z), (30) 

where z = k — v/{2X) and < z < 1. Since the absolute minimum of the function F{z) from ([30)1 corresponds to the 
point z — 1/2, each branch Ick of the critical curve lies to the right of the vertical line A = Ac (in the plane (A, i^)), 
where 21n(4Ac) — i^(l/2), i.e. Ac = sa 1.78. All the branches of the critical curve Ic as well as the phase portrait 
of the initial model in terms of (A, u) are presented in Fig. 1 (left picture). Clearly, there is a strip < A < Ac which 
lies, as a whole, inside the region, corresponding to the pion condensed phase. 

In the right picture of Fig. 1 the behavior of the pion condensation gap J vs A is depicted aX v — l.h. It is 
easily seen that this quantity oscillates as a function of A. However, the amplitude of this oscillations is a rapidly 
decreasing function of A when A 0. Similar oscillations of different physical quantities such as gaps, critical curves, 
particle densities etc vs A were also observed in some NJL-type models with one compactified space coordinate, but 
in a qualitatively alternative case with nonzero baryonic chemical potential [30| . Moreover, oscillating phenomena 
as functions of curvature are inherent to NJL models in the Einstein universe, i.e. in the curved spacetime of the 
form X [3l|. In Fig. 2 the behavior of the gap (5 vs is depicted at A = 1 (left picture) and A = 1.7 (right 
picture). Concerning this type of oscillations of the gap (5, it is necessary to note first of all that its period is equal to 
2A. Moreover, it is clear from Fig. 2, and this fact is supported by numerical calculations, that the amplitude of the 
oscillations of the quantity 5 vs is a very slowly decreasing function of v. Finally, it is evident that the smaller A, 
the smaller the amplitude of this z/-oscillations of the gap 5. 



B. The case of antiperiodic boundary conditions 

In this case a = 1, so in instead of Vl one should use the effective potential (|A12p . Then in terms of the 
quantities (|24p we have 



Da.(to,5) = (m^ + ^2)[ln(A) - 7] - 2A ^ | ^ (^./^^^{2n + 1)2A2 + + S"^ 



+ \j {^^Jm^^ + (2n + 1)2A2 - + 8"^ - 2{2n + 1)A 



m2 + 52 
(2n + 1)A 



(31) 



The critical curve Ic which divides the parameter plane (A, v) into a region with symmetric phase and the region, 
corresponding to a pion condensed phase, is now defined by the following equation 



8(6^) 



°° f A A 2 1 

P=ln(A)-7-^| ^^^^ ^ ^^^+ |/2j^+i)A-j.| " 2^rTTj " °' 

n— ^ 



(32) 
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FIG. 3: The antiperiodic case: The phase portrait (left picture) and the gap S at v — 7.5 (right picture). The notations are 
the same as in Fig. 1. 



As in the case with periodic boundary conditions, for solving the equation p2p it is convenient to represent the 

oo 

parameter (A, i^)-plane as the union of Wfe regions, {X,v) = [J u>k, where 

fe=o 



wo^{{X,iy) ■.0<iy<X}, ujk = {{X, i^) : {2k - 1)X < < {2k + for k>l. 



(33) 



Accordingly, in this case the critical curve Ic is composed of different pieces, i.e. Ic ~ [J Ick, where Ick is the part of 



k=0 



Ic, arranged in the corresponding region LOk- Obviously, the equation for Ico is just the equation (|32p with omitted 
absolute value symbols. However, the equations for Ick {k > 1) look like 



Irk 



fe-i ^ 
ln(A)-7-E 



A 



A 



{2n+l)X + iy iy-{2n+l)X 2n + I 
~~ A A 



-E 



(2n+l)A + i^ (2n+l)A-i^ (2n+l) 



(34) 



Summing in p2p or p4p with the help of a program of analytical calculations, it is possible to obtain the more concise 
form of the equations for different pieces Ick {k > 0) of the critical curve: 



Ic 



2 ln(4A) + 2^P[k 



V 

2X 



2X 



+ 4' 



V 

2X 



= 0. 



As in section fill Al the equation (|35p can be transformed to a formally Wfc-independent expression 



hk 



21n(4A) = -2V'(z) - 7rcot(7rz), 



(35) 



(36) 



where 2; = fc + 1 — v/{2X) and 0<z< iforfc = 0, whereas < z < 1 for /c > 1. Note that the equation ([M]) coincides 
with ((30)) except for the different X dependence of the variable z. In Fig. 3 (left picture) the first several branches 
Ick of the whole critical curve Zc, which divides the (A, v) plane into a region with pion condensed phase (the number 
2 in the figure) and a region corresponding to a symmetric phase (the number 1 in the figure) , are represented. Note 
that the strip < A < Ac of the plane belongs to the region 2 with the pion condensation phase. 

In the right picture of Fig. 3 as well as in Fig. 4 the oscillating behavior of the pion condensation gap (5 vs A and, 
correspondingly, vs v is depicted. The properties of these oscillations are the same as in the periodic case. Namely, 
at fixed v the gap 5 is a quickly damping oscillating function of A when A — )■ 0, whereas at fixed A the gap 5 oscillates 
at — > 00 with a very weak damping. 



IV. SUMMARY AND DISCUSSION 



In the present paper we have studied the phase structure of a two-dimensional GN model at nonzero isospin chemical 
potential /i/ and in the spacetime x with nontrivial topology, when the space coordinate is compactificd into 
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„J . . . . , . . . . , . . . . , . . . . , . . ^1/ .J . . . . , . . . . , . . . . , . . . . , . . . . , rl/ 

" 2468 "2468 10 

FIG. 4: The antiperiodic case: The behavior of the gap (5 vs at A = 1 (left picture) and A — 1.7 (right picture). 

a circumference of a finite length L. The consideration is performed in the leading order of the large- A'^c expansion 
technique. ^ 

It turns out that in the case with L — oo the pion condensed phase is realized in the model at arbitrary nonzero 
values of /x/. In this phase the corresponding order parameter, the pion condensate A, does not depend on the isospin 
chemical potential /i/ and is equal to Mq, i.e. to the dynamical quark mass in the vacuum. The same phase structure 
at /i/ ^ occurs in the chirally symmetric QCD, where pions are massless particles, so one more common property 
is found which is shared both by the GN model and QCD. As a result, the assurance that the finite size {L ^ oo) 
effects of the GN model are also inherent to compactified QCD at /x/ 7^ is raised. 

If L is finite, then the phase portraits of the model in terms of A ^ 1/L and v ^ [ii are found for the case of periodic 
(see Fig.l) and antiperiodic (see Fig. 3) boundary conditions. Among the most interesting properties of these phase 
diagrams is the fact that the strip < A < Ac « 1.78 lies as a whole inside the pion condensed phase. We have shown 
also that the pion condensed gap 8 is an oscillating function vs both A (at fixed v) and v (at fixed A). The same is 
true for other thermodynamic quantities of the model such as pressure, particle densities etc, and is inherent also to 
the (3+l)-dimensional NJL models with curved spacetimes [3lJ or spacetimes with non-trivial topology [soj . 

One more interesting issue should also be mentioned. Although in the present paper the spatially uniform pion 
condensation was assumed for simplicity, however, for a sufficiently high values of /i/, the pion superfiuidity with 
inhomogeneous condensate might be realized in isotopically asymmetric and spatially infinite quark matter systems 
(3^ . A detailed investigation of this possibility in the case of finite space volume is outside the scope of this paper 
and should be left for further studies in the framework of different QCD-like models including the GN model. 
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APPENDIX A: EFFECTIVE POTENTIAL IN THE VACUUM (/xj = 0) 

Note that in the vacuum and in the i?^ x B} spacetime, the expression Vb for the effective potential of the initial 
model (H)) can be found starting from the TDP pT|) at /i = /i/ = 0, where without loss of generality it is possible to 
put TTa = (a = 1, 2, 3): 



(27r)2 



In 



/ 2 2 2\2 
(Po - Pi - cr ) 



dpi 
27r 



Pi- 



rn 



^ It should be noted that the problem of the IR-behavior of the correlation function of quantum fluctuations in two-dimensional QFT 
models was extensively discussed in literature with relation to the Coleman-Mermin- Wagner theorem. One may mention the papers on 
the 2-dimensional GN model [Til . |15| . [T^ . ITtII , where this problem has been investigated and it was demonstrated that for the limit of 
infinite Nc this theorem is not valid, and hence spontaneous symmetry breaking may take place. 
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Here the second equaUty is obtained with the help of formula (fT3)) . Since in (jAip the last integral is an UV divergent 
one, we regularize it by cutting off the integration region, i.e. supposing that \pi \ < A. The effective potential 
Vo((t), as a whole, must be a finite quantity at A — > oo. So the UV divergence of the integral term in (jAip must be 
compensated by the term with the bare coupling constant G, which, of course, has to be a A-dependent quantity. To 
find an appropriate expression for G, let us recall that in the vacuum the chiral symmetry is necessarily broken down 
in the framework of the model ([T]), and quarks acquire a nonzero dynamical mass Mq which is a nontrivial solution 
of the gap equation Vq (cr) = 0. Taking into account this circumstance, one can immediately obtain from the gap 
equation the expression for the bare coupling G. Substituting it in (|Aip . it is possible to get for A >> Mq the 
effective potential ((20|) of the initial model in the vacuum and in the x spacetime. (More details about the 
above renormalization procedure for Vo(cr) are presented, e.g., in j2ll. [23j). 

Now let us find the effective potential Vl (c) of the model ([T]) , when the spacetime has a nontrivial topology of the 
form R^ X and quark fields obey the most general boundary conditions (PT|) . In this case one can start from the 
equation (jAip . in which it is necessary to perform the euclidian rotation [po — )■ ipo) and then use the transformations 
according to the rule ((22|) . As a result, we have 



oc 



n— — oo 



dpo 
27r 



In 



(A2) 



L 
2^ 



cos(n7ra) 



Let us next use in (|A2[) the formula Ina = — — e which is valid up to an infinite constant independent of the 



parameter a, as well as the Poisson sum formula [2£ 

oo J- j — oo 

n— — oo n— — oo 

After integration over po, one can easily find 

1 °° r- 



ds 



cos(n7ra). 



where 



^"('^^ = 4G + 2^ I 



(A3) 



(A4) 



(A5) 



is another, equivalent, expression for the effective potential (IA1|) . (It is clear from general considerations that at 
L ^ oo the effective potential Vl(ct) must coincide with the effective potential in the R^ x R^ spacetime, i.e. with 
Vo{a) given in (|A1|) . Looking at the formula (|A4p (or (|A7[) below) at L ^ cx), one can easily obtain in the right hand 
side the expression appearing in the right hand side of formula (jASp . So, due to the above reason, it is the function 
Vb((T).) Moreover, in the following we will use for Vb((T) its renormalized expression ((20|l . Let us integrate in (jA4p 
over s using the well known relation [32i | 



A 



n/2 



= 2 ( - ) K^i2VAB), 



(A6) 



where Kn{z) = K^n{z) is the Macdonald function [29]. Then, 



TT , ^ Tw N 4(7 ^ cos(n7ra) , ^ \ 

IT ^ ^ " ' 



TT — ' n L 

ri=l 



Due to the relation z-j^Ki{z) = —Ki{z) — zKo{z), it follows from (|A7p that 

Vl(ct) = — Vo(cr) E cos{mTa)Ko{nLa). 



da 



(A7) 



(A8) 



71=1 



The series in (jASp can be modified appropriately [s^, so (here we use the effective potential Vo(cr) in its renormalized 
form (HOD) 



d_ 
da 



Vl{ct) 



2a^^/MoL\ 2<T7 



2a 

TT \^ 47r y TT ^/oHX+n^a^ 
1 1 



2-E 



^7r2(2n + a)2 + L2^2 2mT 



2a J: 



^7r2(2n-a)2+L2o.2 2n7r 



(A9) 
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where 7 = 0.577... is the Euler constant. Integrating both sides of (|A9p over a, we obtain the final expression for the 
effective potential Vi,(cr) of the initial GN model in the vacuum, when the space coordinate is compactified 



2 

1? 



- —rVcr^L^ + 7r2a2 + — — 



E 



2 r2-l 



V7r2(2n + a)2 +L2^2 ^7r2(2n - a)2 + L2^2 



4n7r 



2ri7r 



(AlO) 



In spite of the fact that (|A10[) is valid for arbitrary < a < 2, wc find it more convenient to have another (equivalent) 
expression for the function Vl{u) at a = 1. In this case one can start again from the relation (|A9p . in which in the 
second series it is necessary to shift the summation index, i.e. n n+ 1. Then, manipulating with convergent infinite 
sums, we obtain 



TT ^- ' 



2ct (MqL 

TT \ TT J TT 



71=0 



1 



1 



v/7r2(2n + 1)2 + ^2^2 (2n + l)7r 



(All) 



Now, in order to get the effective potential, one should integrate both sides of this relation over a: 



TT V MqL 



2 

CT 7 



L2 



E 

ri=0 



a2L2 



V7r2(2n +1)2 + 2.2^2 _ (2„ + 1)^ 



2(2r 



1)tt 



(A12) 



Finally, a few words about the phase structure of the model at finite L. It is clear that if the antiperiodic boundary 
conditions are imposed on the quark fields, i.e. a = 1 in (|21[) . the parameter L plays the role of the inverse temperature 
in the ordinary GN model, but only in the vacuum (with zero chemical potentials). In the last case the critical 
properties of the GN model are well understood (see, e.g., in [11]), so, by analogy, we can conclude that at L < Lc — 
jgj^e""'' the global minimum of the effective potential (jA12|) lies at the point ct = 0. In this case the chiral symmetry 
of the initial GN ^ model is not broken. In contrast, at L > Lc the effective potential (jA12p has a nontrivial global 
minimum point. As a result, chiral symmetry is spontaneously broken down at sufficiently high L. 

The situation is however quite different at periodic boundary conditions, i.e. if a = in pip . In this case the global 
minimum point of the corresponding effective potential (jAlOl) lies outside the point tr = for all L 7^ 0. Indeed, it is 
clear from (jA9p that the derivative of this function is negative at sufficiently small values of a, so there is always a 
local maximum of the effective potential at the point (7 = 0, and the chirally broken phase is realized in the model 
for arbitrary L ^ 0. 
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